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ABSTRACT 


A new adaptive finite element technique suitable for 
complex heat transfer and fluid flow problems has been 
developed and its edge over other existing adaptive grid 
techniques has been highlighted. The novel technique provides 
an automatic procedure of remeshing the finite elocnent grid in 
such a way that high gradient regions receive due importance, 
which in turn results in an accurate solution of the physical 
problem considered. 

The grid generation procedure is carried out in terms of 
body-fitted coordinates ; and y , A poisson grid generation 
system is used for obtaining the ? - y coordinate lines and 
the spacings between these lines in each direction are 
controlled with the help of suitable control functions. The 
gradient of the fastest varying solution variable have been 
used as the control function, i*ath the sign of the second 
derivative governing the direction of grid movement. 

The technique has been tested on three problems, namely, 
two dimensional steady heat conduct! cai with non-uniform heat 
generation, reacting flow inside a rectangular channel and 
boundary layer flow over a flat plate with heat transfer. 

The results indicate that the present method works successfully 
for these types of complex fluid flow and heat transfer problems 
and that it is very effective in terms of accuracy and computa- 
tional time as compared to unadapted grids. 




Sine® til® advent aid digital computer nuaaneai 
modelling of difficult engineering prcfelems has beeom a 
viable alternative to expensive experimenial aodelling. 

Recent y&B.m iiave vdinessed a 4tTm<mdm$pLncTeam in 
the speed of computation, wl1ii the result that problems 
i^ose solutiois wero thought to be cumbersome and 
unsolvable hlther-to are now mm solved efficiently* Sevoral 
numerical ^ehniques have been developed of vihieh the finite 
difference and the finite element techniques have gained 'dd# 
usage for finding the numerical sol»rtl©ns of engineering 
problems* The ftnlte element method vMeh provides the 
capability to develop very ^sineisal and powerful pacteges is 
being accepted increasingly, in all dlsCipHiies Including fluid 
mechanics and heat tmnsfer. Both ths 3«»th©<te rec^lre the 
diirlsi€»l of the i:4iysical domain of the problem into a mesh 
i«lth a finite nw^r of intersecting lisies itiich forms the 
coBcmtational grid* The present work is a complement to 
the tbos^ mentioned methods for finding the solutions of the 
transimait nogHliitetr and oti^r ccwolex problems. The work 



It ha# bBm noted by several ati thor# C B J that when 
the grid is ri#t laost miiaerical solution umthod# work well* 

The mjcrortcal errors are reduced a# the grid adapt# to region# 
of large soluticsi variation* Criteria such as orthogimalityt 
smoothness and concentrati cm of grid lines can be used few 
grid adapts ti on in order to get the most desdrabl® grid 
structure* The adaptative 'irld is lost effective #ien it is 
dynamically couoled i^lth the physical solution (say temperature, 
velocity ate,) so the solution and grid are solved for 

to‘?ether in a single' continuous oroblea* 

The various approaches Fcft grid-adaptati cm that have 
been develof^d go far ate essenU*^lly variational methods for 
the extremisEiti cm ei som soluUai p:(porty* Sciveral consider-* 
atlois are involved bore, some of '•/Mch are sasatimss conflicting* 
The points must conisntrate and no region can foe allowed to 
hmom devoid of ooints* The cEstribution also smist retain a 
sufficient degree of sisoothmssi and the gpgid must not foecerae too 
afeei'»d, else the truncatlcii error will foe increased* This meafiw 
that point# must not, move Independently, but rather eadi 
point must sosaa hw/ be coupled atleast to tt# neighbor#, AliO 
the grid point# inist not awe too far or too fait* else o««illa» 
tion# aay occur* Finally solution error* or otbir drlsdn® 
measure* must foe sensed* and there must foe a aethani*® fair 
tramlititig thoffe into tl^ m&tlm of tlie grid* It thould fo# 
noted tlwt the use of in’ adaptive gri d smy not neoesi^fily 



increa#® the ©oispaiteir tiias even though more c<mp«tatifflns 


are necessary inelucJing those for d^tal-iing the grids. 

However convergence properties of the solution may foe 
improved and certainly fev»r points will be reciulred, 

PBEVIOUS ¥imKSs 

A number of stodies on adaptive grids are based on 
the point that the error can be radioed fey distributing tlit 
grid points so that esm positive weight function is 
equally distributed over the field. The varlatlcnal methods 
or the other methods applied earlier are based on finite 
difference methods and selected different expressions for 
the weight fumtlm \%x) which essentially contains 
first ^rlvatives of the solutiai vAth respect to ph^ical 
space, IWooever* tJ^ variational methods t^tain the adapted 
grid, by solving the Euler^legrange e^^ations obtained by 
mlniiaicattcft of the integrals 1^# and 1^ which represent 
paoothnesst orthogonality and etmcentraticn of grid points 
respectively* 

The adaptive grid peoc entire for the finite differwice 
Bie-Wiod has been polneered by Dwyer and SandersClH^ ffieickblll 
and StlJ^sinD^l^ Kai and ^dterscnCsI], i3neff©C4ll and l^lilailMtshl* 
®t 111* C ^ 1 • Dwyer laid Sanders have applied iyb# smcpthn^se 
approath to a iwoN-dl«iensi«nal cuisvilinear coordinate systeai 
with the grid pcdnts c<nstraiiied to move along me family of 
filed cow#r#tia'^ lime in pcoporticn to the grtdd.®nt of 
dependent 



where ^ stands for tite coapiitational spacer s for the i:feysical 
spaoe and T for the dependent varlaolo of the original problem* 

The used or difficult y with the method ic the lack of control 
of ^id skewness. The work by Qmffo is based on spring 
analogy where springs c<»inect adjacent mesh points and spring 
constants are a functim of the gradient of some property afid 
betweMem the grid points. This work Is extended by NakahasWL et.al. 
v.tjO ctmsidered a system of tensiori and torsicanal springs for 
nwiiti<-diraensior!ai grids I'iiich ccntrcl the npacings and skewnest 
of the grid lines, Thp major draw back in the abewe aentiwied 
methods is the finite difference basis m wMch they depend. 

This reepires cew^lax transformations of the eq^jiations fresa 
physical co'«^dinates to computational c e-ordinates and the 
numerical calculations of the derivatives is inaccurate. Alt® 
th« maaerical code may not be straight forward and certain 
parameters in the codte are left to trial and error as seen 
in the «K»rk ef OasDsl. Rai and Andersen proposed an attractitii 
and r*palsi<Ri criterion for grid point movement. This is 
aocosaplithed by aesAgsiing to each grid an attraction 
proporttcsaal to the difference between the mapiitude of tosie 
mtaeure of error (or soluttai vartari«i) and the aveitagi 

I 

magnitude of this measure over ail points. The use of highiir 
dtrlvatlvoe in the error expressions t*as foimd to result ia 
noisy snetturts and erratic point motion* This fafotedure doea 
not exertist any control over either the aacothnecs or er^ofcnallj 
and ^ grid# and thttefore As tortioii is possible* Th^ i« i 



another apporoach, ihe chaialcal reactic« analogy* ^ich hat 
not be®!i studied thorou#ily« 

Thar© ar« recent wrks v^ich couple finite elements 
and adaptive grids. The moving finite elements of Miller 
is a dynamically-aciaptive finite element grid method in viilch 
the grid point locations are made addtitmal dependent 
variables iji a Galerlin f ormulatl tjo * Th? soXiillcai is expanded 
,n piocetf^lsa liniar rtuK^tiafis* in terms of its values at the 
grid points and of the grid ooint iocations on each 

el<^aent« The c sid.al is "hen c^gjirod to foe orthogonal to all 
the basis finicticn^ for both the so'^ntlon and tb& grid* The 
grid point locations arc thus obtained as part of the finit» 
olament solut.lcru An tntemodal vi'-rortty is intro diced to 
pcnaliae the relative isoiion bet*, /eon the grid ooints and an 
intern oc^i repulsli^*’ force is introduced to maintain a f^imiaa 
point separation, both ©f tt^ae effects are strong but of 
sfjost range. A small long ranp^ attractive force Is also 
introduced to kvTO the no<fes more oqtially spaced in the 
absence of solution gradients. fHher recent work# like that 
cf Berger and JamescaiCT/Il refino ihe existing grid at the 
locatiiais where there is solutioi'i ©rrc®. These are also autowsatii 
adaptatiORs but the drawbacks of most of the#® axe usage of 
ccaaplicated numerical c«Kle and an increase in the coraiiutaiicisil 
itwc (bte to the increase in the ambeit of nodsa. 



Th® present work Includes a ncs^el aii^oacli 
uses finite elements and adapts t»ic grid by solving an 
elliptic system of qrjud generation oguatiaiSi, The pr<»©~* 
dure can 0 © applied to <inet two ana three cJiiaensions \id.th 
mm ease* The method is tested on the followiiig three 
problems* 

(i) Tv#o<*climns-’.onal steady state heet conductieii 

(ii) A stead/ chi?rfiically reactisn§' in a rectangular 
charnel 

Ciii) The stead^^ boundary la'i'cr over a flat nlate T-dth 
iteat transfer. 

In all the cases adaptation Is applied on a coarse 
grid and the adapted grid results re compared 7ith those of 
a coarse uniforri gric and a fine onifoin grid. Comparisons 
have also been made vdth exact anrlytical result® where they 
are availitole* Thus ttw^ study h lags fortii the general 
applicability and the suiiabilix’i of using the no^^&l adaptive 
finite eleiasnt for the solulicn of complex heat 

triarasfer, fluid i low and co»i>u&ll -’n /:oblc«s* 



CHAPim 2 

FCB>mATX€H OF THE AI3APTIVE Fi mt B ELHI^.OT mOCWim 


SASIC CCNCSPTS If*! 'd!iflEriATX(»IS 

Grid f^fieraticHa techniques ywrovlda a #|«teaiatic 
procedure for ill© placement of points v/iiMn the solution 
domain so that a nuuierical grid v4th the desired control 
over grid spacing, orthogonality of cp:id linear and 
swoo’Uiness of grid line variation can be obtained, Sc»iiv? 
of these techniques exercise control over the node mwement, 
vitiilst others genf,rate gild lines tisinq compirtaticml coca?«* 
clnates ^ and vihos® intersection points provide the 
reoMired grid. The attracticsn/ re pulsion method of grid 
adaptatlc*! C. ^ H transfiriite Interpolation 

method C 8 3 vshich Involves an intsTpolation for the Interitnf 
nodal pdUits usinc the preserlfeed coordinate informatitm for 
the boundary n^es, both result in algebraic grid genetatiiwi 
equatisns* On the other hand# the methods vslitah generate 
tl» aoii^tatioiial eoordinaie lines 5 “^onst aiid'^-wcoastj 
invaHably use partial differential efl^ti<Mis af ellAptlt# 
pai^bolie or hyperbolio ty^# Tlwy ai« to find the apinfopxlatt 
transformation from the physical coordiaates x#y ( for a 
problon) i^e compw'feational coordtaates ^ aiKl The 
solution of the i^ysital problem itself way be .carried «Mt 


«ithdr in of tt» ptiyslcal caordir»ates or the corapti*^ 

tatimal coorcJlfiataa* vMch evor is convanient. 

In the present viork^ the ellip-tic grid gsneratitm 
sVfetaa^^fihich is one of the most popular techniques^ has been 
employed in conjuncxim \tlth the finite eleraeni method. If 
FEM is used for solvinq the piiyslcel problem also# ihm 
the coaptitati«ans can be conveniently performed in textns of 
the physical coordinates theiaselvc'S* instead of transforming 
them in teims of the compitationeX coordinates*. 

ELUPTic sysra?? 

The adaptetion of th^' grid can be ¥i&md as the 
mwement of lines/nodes due to the solution gradients .ie. 
the lines/nodes are moved towards hi<fh gradient regime and 
concentrated* Here use has been made of the properties 
of elliptic grid generatim sysWui* Soias elliptic systems 
exhibit the extranma principle st i*0* extremsusf soJutimi 
cannot otewr uithin the field* Such a fea^tore is useful 
for grid genetaiion since tuen maiCiiam end skininum coordinates 
oi a doaaifi occur m. the bounefery* Another liaportmt 
property in regard to coordinato system generattoti is the 
inherent tmooiiiness that prevails in ^ tolutim of elMptit 
systems* Furtlwriaore^ boundary slope i^scmtiiiiiitles are 
not propagated ihto the field* The smoothntts pri^rty and 
the extreraw p«inciple«t allows grids to be fsneratted for mr 
conUgoratiCii without overlap of grid Hues# 


POIsroM SYSTEM 

Control of thB eoordinat© line distribution in 
the field can be exercised by generalising the elliptic 
generating system to a set of poisson equationsi 

(i» 1,2,3) (r>4) 

in vA^eh control fmctions p caa be fashitxied to control 

i 

the spacing and orientatim of tho coordinate lines* Here ^ 
are the eexaptitationai coordinates* Tho extrerau® prl'ntilplas 
may be weakened or lost completul, ylth sudh a systesi, but the 
existence of an extreiuiim principle is a sufficient, but not 
3 necessary, condition for a eonc-co«ofi® mapping, so that 
smm latitude can m taken in the form of the cmtrol 
functlcxis, 

ATTIIACTIOM OF iiODES TO OOORDIMATE Llffi^POIOTS 

Let us ecmsider the grid geometry shown in Fig ,(2,1), 

2 

For the pnlsson grid generatlcxi equation V -r^ »» Q, negative 
values of Q would tend to make the "^^lines move in tte 
directiin of iiicreasii«f • While negative values of p 
in V • p vdll cause ^ **llnes to move in -Mie direction 
of decreasing ^ « These effects are illustrated in Fig* (2#i) 


for an n^liiie bousitfery* 



If ihet points m tlie boundary are free to move on it 
thm ti^ ^ •alines vJouXd move tcRriard Janes lower values 
of 5 the other hand if the boundbry nodes are fixed 

only the inter seeti angle will change at the boundary* 

Thus om may ccficlude that* in gaioral* the negative values of 
the laplacians of the dep<'nd«at coordinates mnm& the lines m 
vdtich that coordinate is c<»istant, to move in the direct! cn 
in wftiich that coordinate decreases and vice versa* 

ADAPTATIOI CRlTERlcm 

As mentioned earlier the control functions can be mani- 
pulated so that the grid distribution can be done as desired* 
Here the first derivatives of the fastest var'ytng soluti<m 
variable are used as the ccHfttrol ftinctions for adaptation* 
Therefore* the control functions and Q can foe written as 

^ dX 

( ?* 2 ) 

i'; Oj 

ay 

wt^re X iHid y are the i^ysical coordinate# and f is tJie 
variable as the basis for adaptation* Ih the developmwst 
of the present grid adapts tiaa procedure# certain idwis fro®® 
the techniqjie of attraction of pc^nis/lifies using control 
fuiicticils have been ineorperated* The high gf^dlent regioii 
in the erl^nal grid tan be teen at a irmi te «**ich lines or 


points ^owld ©onverg«f* Thia is decided fey the magnitMdt 
of the control, function and th© siqn of v^ile tl» 
magnitude gives the relative movemontt the sign gives 
the directicai to which the points/ 3i.»s ^c»ild move* 
the saagndtude of the gradient is increasing In the increasing 
x-dlrectlcati the cmitrol functions should have -fve sign 
and Vic® versa* 

The sign of the control function therefore* can fee 
derived by constderiSH? tte sign of the secorid derlva'^lve 
of the basis fUnctioni, Here vre mn discuss cases in 
general • 

(i) first derivative -hre and its magnitude livsroasts 

(ii) first derivative -hre snd its s^gnituds ctecreases 

(ill) first derivative ««v@ and Its magnitude increases 

(iv) first derivative -*ve and its magnitude decrease* 

(v) constant first dertvative* 

In all the above cases lencept the last mwi) 
m can observe that 



controls both wignitiide as well as the directiefi of node 
»«nr«i»iit in the VMllrectioR* Sams Is the Caso for tins 
©th«r <at^ti®» y# in a tw«*dli«jaionia caso* Finally m 
odd a «fei#it •W* to t^trel fisictioii for restricting tSu 
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Here ’W* is a -fv^e niiraber ^i’iich may be large. The 
roacon for taking large * ;/• is that othervjise the control 
function on the rif^t hand si^ uBy v»0ak®n or ccHuplately 
vxolate the extreimiia principle of the elliptic system. In 
order to avoid the pQ^etllcameni of grid points crossing the 
boundaries# the value of t\p should be as large as possible. 
The value of however can be decided by the maxirawm 
gradient occurring in tho domain for the solution variable 
used as the basis for adaptation. Another functi^xi of 
la to control the smoothness of grid. The large value of »W» 
will not allm very much coneentrated md rarlfied toma 
in the grid, 

FCEMtJLAtICM FCm THE FE PaWGIs 

The present schem# is based «xi a dT^waiioally adapttv® 
grid, S§0 the grid sh^ld h© rearranged according to the 
variati®* of tJ^ solutiati in the i^yslcal irrdfele## say 
teaperatWPt in a heat trwnsf er ^obiem* mmt ftmU 
dimrm&rn or Unite elements can b# i^#d along with the grid 



a<SaptatXai proeedfeire for solvino tf)ft physical problaia* 

Here finite elsfiente is preferred since it is versatile* 
Moreover it can be directlY applied in the physical 
region pnlll!» the firdte differences /^here transforasaticai# 
xo a rect^gular eoraputatl mal domain are a imjsi . 

The only disadvantage ul‘ih the sche®© is that tvo 
eqitatic»is are io be solved for new coordiaates* But 
this can be overlooted vhen m consider the solution 
accuiady and convergence (especially for norjh-linear faroblesns 
Also the computer time necessary vox- adaptation on a 'jroerse 
grid my not be larger than that for a fine grid of same 
accuracy* This point has been repeatedly mentiimed in the 
literature CS 3 * In the grid genoratim schenis the poissm 
system given by equations (5^*3) Is solved by the Galerhin 
technique or the v^eat forsmlatl on, Thu? we havet 



** 0 ( 2,4) 

where the i^discript n dwetes the value at a particular 
node », The equatiase (2,4) can Jae inte^ted by part® t® 
l^t. 




XI- -SI. 


™ II / (^) dn. ( -^.s) 

TbG line intograls arisimj the integration by 
parts heeorae lero '^hen tlie boyntlar’f conditiCffi is DilUCII i,ST 
t'/pe* The boundary con (fit leas for grid g€Mrt@r»tion are such 
that either tha boundary ccaifornin'j coordlJiate is specified 
or the nomal derivative to the bo'UKiary is siaclt equal to 
Eero* This allo'.vs the grid point’: \g ncfve <»» betindarjies 
also* So* ifeiifes© boundary cciKfiti'^r.a (.rs straight forAtard 
and can be easily incorporatec? ini e tlise grid 
scheiue* 

C3I1 right hand side of ec^iationa C a 5 ) tlw fiipt 
and second derivatives of the solntion variable (#iich is 
used for adaptation) appear and these* oontrlbntlofis need to 
be oaic^ilated at each node* In the evnlttatlon of the eleiwntal 
Integrals# however# it it cetteisary to ym local coerdSjiatet ^ 
and V liiitii fatilitat® nwaerlctl integratic!ii* Mr Ineegnlar 



shaped element, for instance* is locally tran.sfor!a 0 d into 
a s<|uare with ^ and V as the local coordinates. All 
derivatives with respect to the ctlobal coordinates x and y 
need to he trantforwd in te-^'nis of tte local coordinates ^ 
ami V and these t ra ns formatitms aro performed according to 
the relations given by Thompson CBU . Therefore th® 
control functlcxis can be expressed in terras of ^ and as s 


fjj ** (y f - y f )/ g 
** (•* ^ t- X f)/g 

^xx * ^ ^ y y ^ + /■ f )/g 

+ (y^ y • 2y y y 4 y Kx f • x f ) 

+ X *«• 2y yx + x }(yf**yf} 

tyy )/g 

+ (x y'<»»2xxy4-x y)(xf«»xf) 
+ X ^ 2x X X t X )Cy f ♦ y f)) / 

( 2 , 6 ) 



The terms such as x. , ^etc, in the above relations 

b 7 

are evaluated by inaning the transfoimatim between x,y 
( glctel coordinates} and ^ ,1) (iocnl coordinates)* For 
isoparaTOtric elements v^om boundaries are expressed by 
the same order of cirves as tb'^ order of interpolation within 
the element, it is oossihJe to ’vrite t 


and 


X » I. 

y w 


Jg 


3K 



( 2 * 7 ) 


where are the interpolation iunctlom fear the solution 
variables* 

After solving the pedes on equations ( ^ So,) and 

(2 5b), the nao^tatiaaai coordinates ^ and are obtained 
at all original grid points. Then ^ ** constant lines and 
const lines are drawn by interpolating between the 
calculated nedal values of S and ^ , The int eSfsecti 
points between the 5 • c^stant «id ^ ^constant lines 
give the new frld points* In order to control the stability 
and cmurexgenee of the evoiuticvi of the grid, suci^ssiire 
underreiaxati on/ overrelaxation technique is also employed* 
Etetaiis of tlie grid programme pactege are giv«i 

in Appendix t# A - 







CHAPTER 3 


TEST PR0BLEM>»1 

nm@il ®f adaptive finite altmenis using a 

polsscm system Ims been tested «m tb?@e problems. Here it 
is intended to #iOw the practicality of the method for heat 
transfer aiwt fluid tswchanScs probleiJ^s and the difficulties 
vMch arise in nulti'-dimoaslonal cases, A FCBTRAH package 
has b®®a developed based m this technique for generating 
a 2-dimeiuiimai numerical grid v&iich may be used to solve 
the heat transfer and fluid flow fnroblems in conjunction with 
thi FBM prcwsedure or any ott^r stiltable procedure. It is 
however, appropriate to ocdnt out here timt if the prc^lam 
is solved by FSM using S-noded rectangular elements, elective 
and straight f orvaard linleing of the main program ( of tte prd>lem) 
to the adaptive grid iMtItagi is possible* 

00M3UCrFI€»l PROffl-EMs 

The first test problem e«sidered is that of steady 
state in a rectangular slab* The egwitlon 

gowrntnf the pr^lea is gtwi by, 

4 . Q»» » g 43^1) 

^ere T is the tejsperature laid is the rat® of heat g«««ra1»d 
per wg4X yoSMse# In two dlmenaicns the sguatioii mm be srittei m 





WE 0 


( 3 , 2 ) 


(K|i> ^ ^ '»"■ 


F©r finding th© finit© ®l©K»nt sointloi ©f th@ ab©f® 
©guatiofii, Gal®rkin*s weighted sfesiduai approach is used, 
as follov^? 


m 

X 


i 
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dy 


f |f i wi ^ 

' ^ s y' 


dx dy 



<hc dy « 0 


(3,3) 


v^re axe ttm interpolating functions for the teaiperature* 
Integrating th© above equatim fey parts md rearranging it® 
terms we get. 



-^x 




dx dy 


jOl 



-a. 


Q« « dx dy + ^ ^ 

(3,4) 


wt»r# t3be last tern indieates th® lin® integral m@r the feeawdary 
(fi) of the solution do«»i« (0-) and is evaluated at tli@ feoinidario* 
for Heiaaiia ©r Mixed feouiidary ecsidi lions. Expanding the yariafele 
T intena# of sliap© funotions interpolation l^ncticns) the 
eg^tion l»e@on»« 
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m 

In matrix form* fco: c«ie elemenrfc the left hand side and 
right-hand side *5onxribyti(K!S( j^esnectively 



*12 


*18 


•“tl 

®21 

# 

*22 

» « a 

*28 


^2 

• > 

m 

i 






; 

*81 

*82 

# e * 

*8F 





For all the elomentg similar LHS and RHS oontriboticxia 
arise wliicli ean he assaaiblsd to form a set of ^algebraic eqmtlms 
for all nodal variables* Xhe final soltition can be obtained 
after snbstilwting Oirichlet condition (if there Is any) by 
matrix inversion of the nodal eqyetions. 


WRlFICATIdit 

The progrstt denreioped for eondooticai has been tested 
for eorri<i’fc!iess «dth the help ©f a teonasu solution* Steady 
state hisit siiMhitticii In a reetanpils^' slab is eoB^idered eitb 



the following data«# as shovii in Fig» (3»l)i 



hij* 

1*0 

%£ssO 


0 

^yaO 

* 

0 

^X«1 


0 

T 

V*t 


0 


(3*8) 


irEesults are tabulated for 4 clen?nts and |6 elenmts 
(both S-4ioded rectangular) along ti.th the series solti-Wlcm* 
The 16 elOEient soltition series scluillon ar® almost 
Identical* See Table 3*1, 


AimPTATXC^l 

In Ihe laesent problem (Fig* 3,2) to highlight the 
effects of grid edaptatloiii the heat generation term is tahen 
as an exponential ftnctiaa of co-ctrdlnatee so that high temper^ 
at«r® gradiwtf occur in certain regi«s. It Is given by, 

Q"' « qT + (3.9) 


The boMndary condlticms aret 



^■^*0 

m 0 



* 1*0 



• 

0 


X<K$ 

m h(T *• 

(3*1 0) 


the pripirty yal-ws and other ton stents in ^e e<itiatiiii« ax# 
taliii m feUemi 



T 

T(o,y) 

4 eieMnts 

"J”! Of y) 

M elementfi 

T(o,y) 

saxlea aolotlen 

0.0 

0.2941 

0.2946 

0,2947 

0*25 

0.2790 

0,270-5 

0,2789 

0,50 

0,2292 

0.2293 

0.2293 

0,75 

0.1396 

0,1397 

0,1397 


f *00 


0.0000 


0,0000 


0.0000 




h « 1 0 

* m 

k « 200*0 

P « to (3*11) 

Tl» pac'csblem is soi^'&d for 2b aiesents (96 nodes) and fcflf 
ICO etepiemts ('^4f nsino vtlfom mesh* The adaptation 

is carried oi*t cm a coarse grid (25 olements)* For adapratlm 
the connectivity^ ^ and y co-ordinalus of the nodes and tim 
soiuticm obtained by coarse grid arc reguired. As mentAoned in 
Qiapter 2 the new co-ordinates ^ md. ^ shmaid be boandary 
ccmforming* The ^ values a4:» « f and x ** lines are 
specified. Similarly values on y « 0 and y ** y^j, are 
specified. The pd.nts are free to move on the bouixtarles als®, 

the specif t«i ^ and'^ hmmm tJirichlet botindary ecnditicnis 
for ii» system of grid generaticm equations (Ecfuaticn 2„t)# 

After #a#ii adaptattcm “Wie temperatures in the slab are found by 
the cimciuttiois progress 

Fifwre (3*2) to Figure (3*5) ah®v the evolution of 
the grid fnxa the initial. UBiiforirt mesh to the final adapts 

The regions of high gradients can be clearly observed 
from these flgwe®. The grid adaptation is ptrfor»d four 
tlwss end the final results cf the teaperatoe vaxiatim is 
presenled i» Flpire and Figure C3*t),fer different 

Holies sahieh lie csi a &hmm curve* For e^aparissn the coarse pad 



adaptive grid aad fine grid rewilts for correap«mdipg 
points ar® plotted l»n the stra® graph, after interpolating 
the coarse grid and fins grid solutions* The figures show 
that the adapted grid soiutim falls hetmm the coarse grid 
soiutiiMi and fine grid solutim. The nearness of the solution 
of adapted grid to that of fine grid is apparent after repeated 
adaptations, 

mm CONTROL 

Definitely there should be sicsae cmtrol over the grid 
point moveTOnu so ars 'cq actiere to tlia principles of smoothness 
.-ind orti-jo^jonaliry at the bounda'Cies of the grid, Truely 
much study has net been dgne on this aspect. But in gwierjil 
ti^ hlgflh value of <W* (+ve) in Egn, (2*3) 'jives stiff resistance 
to grid point movement. This is te the reason that hic^ 
value of ’W* m^s the rt#it hand side less and in turn the 
Laplacian projaerties are retained* For initial guess the value 
laay be of the carder of absolute marlnuier. of the solution gredlent* 
The value of is increased fot each successive adaptation 
by a mmlX amount so that siaoothness trouble may not arise* 

An iiipoa?t«nt paint ^hat needs to be mntlmmd is that the new 
grid points say shoot out of the ooundariet if the value of W 
Is very eisall Kdien compared to the gradient values. In the 
current problea the value chosen tor 'M is 9>©*D and it is 
increased t*25 tims aft«.>r each adaptation* Fri» the above 
digcussioit it is cle^r that the parameter it problem depaftdent 
mt4 vaxleiiis acasrding to ^ mximm graditnt )Cai>s« value) 
of solutloii^ 









1 

I 

f 


I 


FIGS 4 GRID AFIER 4TH ADAPTATION 




FIGS. 5 ISOTHERMS FOR ADAPTED 25 ELEMENT GRID 


nn 



F1G3.6 COMPARISION OF TEMPERATURES FOR THE ADAPTED 
AND UNADAPTED GRIDS AT CHOSEN NODES 




FIGS. 7 COMPARISION OF TEMPERATURES FOR THE ADAPTED 
AND UNADAPTED GRIDS AT CHOSEN NODES 




CHAPTER 4 


CCMBUSTIOH PROBLEM 


The second problem considered demonstrating ■Wie 
capabilities of propo^sd adaptive finite element technique 
is the lamiiiar flow viith chemical reactioi! in a rectangular 
channel , Ihic ic a t'^ical ooBhuatim prts^l&m #?ich 
is highly om-X^near '"ith sever© ec«v®i^g®nce difficulties. 

Combustion proolesm eusually e^re characterired by fast 
chemical reaction processes #>ich load to high teT^fjerature 
and concentration gradients in certain regions* the grid for 
nuisaarical modelling rust be suXtah]’/ formed to capture these 
hi^ gradients falling v^iich sevorf' inaccuracies in the flnaii 
solution may occur* The lotatXcst and magnitude of tha high 
tPfnparatore and concentrati c*i gradients are not available 
a priori and therefore, adaptive orid prosedir® is a v«ry 
effective mt to tactic such prdbl©?3S* 

The present jiiyslcal profoler treated as a steady 
two-dimensional flm of gas in v^ich transport of throe 
chemical species - a fuel, an oxidant and a productpoccur. the 
reaction betisjeon tho oxidant and •f’uel rcltases ht#t and this 
hast is also tranaperted bv the fluid flow* Tt« loosl chtnpit 
In temperature «id species concentration modify tlie local 
(dtensity md hence the flow pattem Itself* TTait it Is «ts<^tlal 




to solve the coupled equations of mass* raomentuv'^energy 
and species to obtain velocities* pressure, temperature 
and concentratim®. 

The governing equations are taken as follwst 
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(4*2) 
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( 4 , 3 ) 
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(4,5) 



The ia^ortant aeeuraptione involved In i^taining the 
above ecmaticats are* 

(1) All ?^lrs of species have tacio valwe of diffusion 
coefficient* 

(2) Body fcKTces are abs^t* 

(3) Pressure wot k and vtsccms dls'j'pation are negligible* 

(4) All prcpertles except density a re^ constants* 

As fpr the third special one obtains fxm 



C4.6) 

S, + Sj + s, - 0 

(4.7) 


vAiere tl^ syiibols x, y* «» v I, ? p k * p and stand 
fcr the usual notation and B stands for the binary diffusion 
coefficient* 


Finally closure for the above set of equations reqpims 
the specification of the mass and heat sources from the rate 
of chemical reaction and iim thermodynamic ecpatioiis of states 
illicit relates density# pretsurst tupecatitre and the mas* 
f^ctiieis* 


frcMi tl# perfect ps law for ideal adrture* 






Hate of prodijetien of species 


I 




-E/RT 


Q V V 


OC 



M, 


) 


</. » 1 


% 


3 

H 


P n 


(4.9) 


(4,f0) 


As fcff «« tl» Pcumidary ccsiditims are cmesisisd or 
solid wall boundaries the velocities have to be set to aero 
but the heat and m$$ ftacilm variables raay have either 
specified flapies or specified values-, On axes of synrasetry 
noriMl velocities and iioryial derive lives of other veri«4ble® 
have to be set to xejco. The inlet profiles of UL * v, T and Xc 
are specified. At the wtlet the velociiyt tenperatore and 
ccncentrstioi profiles see taken to be fully developed. 


THE FIHIXE EUKlEht HOUATICMSt 

The program incorporated isomrauntric tleiaents cmsls- 
fcing of quadrilaterals faving c^adratic variation. 

The pressure is represented by linear functions because of 
the Xmmt mfd»r derivatives of pressure and Idie regaining 
v*rieble«f «# v* t$ T end are represented by quedrstie 
fbnctiemi# The Salerkin technique is applied to ^gcttJsc 
finite eleaent equations* For eraaple the epecies equatim 
14,»5) 0m be ietefrated at feUM-i 




wJier® p the ship© fwneiicns for q'ndratic variation md SL 
i 

is the soiutlcai diMaiii* 

Or th« part of the betmdary the variable is 
specified and m ih& boundary T 2 the derivativt- bY</en is 
specified* 


The other eqwatiais can also be integrated exactly in 
a tlmilar wiy* The contimiity ecpiaticm is bandied differently 
as dhsert^ed in other works in MteratwreD^ 3 • This is integxatec 
with resoect to linear weighting functions as follows! 



dCSh) 

S"T^ 





<ja. US 0 

i w 1#**4 


C4*t2) 


vdiere sha|w fhncticns for linear variatiesi* 


CHAN«. m.m$ 

Two r#iE#tiii® specie# flow along n tncMi inensl ina 1 
thnnnel and hest released disitirhs ^ flow ^«h it 
ijMi Inlet h#®* Iteily develoiied# cssst of tiw species* <4* f imy 


reaction laay be regarded as be lag 0 :fitectlvely of first 
order and species e<Hiitions other than that of need 
not be eonsidtred* 

The source term for the mats fraetion of the diiiiW 
s pedes » appearing in equatioB 4*5 ir then tahen to be 

^ ^ C4*I3) 

^i-^ich it independent of the temperature* For the eguatias 
of state VII® asewmr the form 

8 =-- 30Q/T (4.14) 

eorrespcndiiig to Wiifows chemicaX cafflpjsitton# For the case of no 
heat release (h^»» 0)t the flovii re»’wdns felly developed ai»d 
the equatim for need only be &oZ\od* 

With ^ origin of ^ rectangular eo»ordinat«m! being at 
the midpdnt of the inlet* the boandary conditions for the 
cases of no teat release and teat release are talcm as 
folloiws (Figure 4.1}# 

(1) Inli' (< - Of -• 4 y ^ a) 

* ” *o 

y «. 0 

t m $00 

f(| W X| 



(2) TOP and boittj® wall® (o £ x £ 0*5, y ** i, a) 

ay, 

« « V W ** 0 

T ** 300 


(3) {x » 0*5t -a^ a) 
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The foIl<mlng values ©f are talseni 


® D « 10*^ «% 


K « 0*1 W/i%?k 

Cj^ » 1000 J/Iio/K 

hj • 5x10^ J/Wil 

h^ w h^ w 0 

«• 1,8x10*^ feg (ai/®) 

!X #■ 1 


ft «* o«05 a 
Lftfigth m 0,5 m 
Width 0,1 » 
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REI^TSte 

gytdL 

te the ftfow neatiaafinl pareblftm it 
{Stained in «r4it to Qe£ied%aft«t ^ th* prei^ttn 

««d« «f the ittidf by ne^ariag lAtti wailAblft 

i^C ^ lidio haft iit»® W.i pfdHai ad^g a 



miFORM ffilD SaLlJTIOMSS 

The do«i„ is divided into nnifors, 8-neded reetengulsr 
oloments with 17 elements in x-directlon end T-ejements 
in ynilreetlor. The total number of nodes is 406 anCtthe 
total degree, of fteedcm ,t, 1768. For numeiloal integration 
9-polnt gauaslan niodretute ha, been used, pieard', iteration 
nyo««dir« is used 7 tiws to »oo«re coovergeneo* 

In th® «#s@ of flow \dth m hast release the csswergence 
is o^t.■=^lnea in the first iteratim iteelf since U., Vi and T 
-rc constant. The ^ofiles cf wss fnetioR at different 
sections has been px -sen led in Plq, (2,a), These caai be 
compared with the profiJes in L ^ which have been 

reproduced in Fig^(2,b^) 

in the case of flow with ha«t release the 
at differmt eections is plotted in Figisxe C4,3a). These can 
be compared '^e profiles giv«fi inL H 1 which have been 
reproduced in Fig, (4,^), The profiles ©f mess fraction md 
velocity shew good coincidence vith that of C ^ 1 * 

The isothexxfe for "dte hest release case are plotted in . 
f^ig*C4*4a)wiiich are coapured tsdtii the isotherm plot ofL^ !§ 
reproduced here in Fig*(4,4b^ 

AQ4PTATI<» 

The tdaptaticsi is i^rriod oat m a saaller grid «n449 
elesitiits Cl <1164} and t4g nodes ^th <»i de^es if 


Gjfadiints ill species heve been cbosen as ihe foisis* fox 
tdaptaticsES vti@n there is iio heat release since the prsfolew 
is Isotiierniil* This cheice for adepfaticsi vd.li smoth 
the discrepancies arisin-i due to the sp^^eies varieticii. 

For the sictiaticsi v^ith heit release a«feiptati«i is carried 
on the basis of tesinerature bteadse of the rapid variatiems 
in toi^peratyre which disturb t!w de^-iity and eeaiseciien tlyt 
the velocities and species. Totally ^iJC iterations are 
carried out for- obialnlng the recfuired result. The profiles 
:>f species and velocities are given In Fig, (4,5) and Fig, (4,6) 
The isothersas are plotted in Fig, (4,t), 

The adapted grides are shovai in Fig, (4,8) and Fig* (4,9) 
The adaptation not only reduces the c<»|3«tati€Bial tiiae but also 
gives the Istoi of the variation® in the fleSd varlablss, A 
renarlable feattire in the idaptation of conbustlQn probleB# it 
that the grid itself shows the location of heat release and 
flaaie. This ««! be i^served in Fig, (4,9), 

oomxtmcMm 

(I) For starting the iterativ® protedusre for ihelilwat 
release case •^e s^wtiaii fro® the no heat release 
east is titon as the initial gaest. This entires 
fast tORvergmee, 







solce of the weiyi-.t *'7* vMch controls the xat@ 

.d rawoTOnt b'^tv^oen t'vo &nceagsivetiiaKiati<»is» 
on ih'* rexiwjffi derivative of the variable nsed 
oasia f'^r edaptatim* Larger the derivative ^ 

^ the value of ^W’ viltich will give a stable 


M&i. solution* 
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FIG4 1 UNIFORM GRID WITH 40 ELEMENTS AND BOUNDARY CONDITION 





FIG4 2A Y1 PROFILE PLOTTED FOR 119 ELEMENT GRID 
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PROFILE FROM REFERENCE [9] 

BROKEN LINES FOR THE CASE OF HEAT RELEASE 


FIG. 4.2B : 



FIG4.3A U-VELOCITY PROFILE PLOTTED FOR 119 ELEMENT GRID 
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FI84.5 Yl PROFILE PLOTTED FOR 40 ELEMENT ADAPTED GRID 



= 0 0623 



FIG4.6 U-VELOCITY PROFILE PLOTTED FOR 40 ELEHENT ADAPTED GRID 






CHAPTgj 5 


Th# lardJi^is solved by the «ppllcatiQH of the 

adaptive finite eleraent teehnicpj® is the Blasins problem ©f 
boundary layer flmi over a flat plate Fig. (5,1). Apart from 
the velocity boundary layert the thernial boundary layer has 
also been included for analysis* 

The j^ysical desenption of the prc^lem is as follows* 
A viscous fluid flcK’s steadliy over a flat plate at zero 
ilielidiiice* The flow is Iggalnsr and there are no pressure 
gradients in the flow directiim* beat transfer occurs due to 
an applied temperatui^i difference between the fluid and th® 
plat®. The equatiORS governing tho flow and th© b&at tensfer 
(for incoiaprestible flow)i 

C5.2) 

(%3) 

vAmm C..# b# o^tstani# and are epeeified for the 

selittien of u# v ewi f m4 f 

tislmf 

r m 

f f/s 
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^ ^ » © 

“i * * ^ ‘ 

30pii||+S<^v|| (t|f) 



and f 3= T/x 

' 00 


( 5 . 4 ) 

th© diaansionlass governing eguationSi, after dropping the bar# 
over the non-dliaetisional qwantities for convenience* beconia 

>* 0 


(5*5) 

It ^ + V • 


(5.6) 

II y m 

bt d y 

JL Si^ 
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WJHERICAL SOtUncm 

Al*lli«»gh a siropla aiffiilarity soiiitlGHi exift# for 
this prel>l«% a finita ei«itnt nwerlcai tojMitlsfi using a 
itnifoxn nash has ais© bmn (^tatnad in ©rdaf t© bring f«wth 
th» advaniagas ©I adaptation* Tha FS5 prcffiedtirs is basad 
cn tha Qalsjrkin wilgh1»d residual approach* The rasldu® 
©cmatitais are obtaifiod as foil^ssi 
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( 5 * 10 ) 


BOIMJARY COJOITKms 

On th^ flat plate #ie velocities are equal to taro 
due tP no slip condition* The tesiperatur© ©r fl«X| 
specified at tlie flat plate (or v«ll), m ^e far stream 
the « velocity beeoittes equal to the free str^m velocity and 
similarly T is equal to At the leading edge u and T 

are free stream value s aad v is equal to fero* ^nce the 
equations are parabolic th» dmm stream emditions need not 
bo specified* These are shcfsei in Fig* (5*1)* 

The retttlts for the wm solution of Urn squatlois 
are plotted In Fig* ( 5,2) al«wtg with the Slasitis series solutlcn 
for cwparisfiR, 

Adaptatisn In tM-s pr^lea Is carried out in two 
different wsyi* Once the grid it moved according to 
the qr^Wonic u «d anolher Um ^e h*^s 4» tomperature^ 

m both the casei the wiifo» grid coJhitlon* adapted 
iolittcd and ilaalue series soiotistt are plotted* im 


(5*2) and Fig* C5i*3) tht Uj-irolocity praPila for the 
u-adapted and T-adaptad solutlcffis art pretented respectively. 
In both the cases the solutim accuracy of *ui? is increased. 
But Fig, (5,2) shows that the u-^dapted soXutim gives 
better result for »u». 




Table 5 >2 


ETA 

u 

unifcarm 

grid 

11 

u-adapted 

grid 

u 

Blasitis 

0.0 

0,0 

0.0 

0,0 

0.4 

0,127729 

0,130326 

0.13277 

0 . s 

0,2557 

0.261670 

0,26477 

1.2 

0,3806 

0.389625 

0.39378 

1.6 

0.5 

0.512254 

0.51676 

2.0 

0.6121 

0.625445 

0,62977 

2,4 

0,711 

0,725592 

0,72899 

2,8 

0,796 

0,8W606 

0,81152 

3,2 

0.864 

0,875982 

0.07609 

3.6 

0,915 

0,92490 

0,92333 

4,0 

0.951 

0,95830 

0,95552 

4,4 

0,975 

0.976 

0.97507 

4.8 

0.988 

0.989 

0,98779 

5,2 

0,989 

0,991 

0.99425 

5,6 

0,993 

0,994 

0,99748 


Table 5«3 


ETA 

u 

imifona 

grid 

u 

T-adapted 

grid 

u 

Blausius 

0.0 

0.0 

0,0 

0,0 

0.4 

0,127729 

0.12864 

0.13277 

0.8 

0.2557 

0.2577 

0.26477 

1.2 

0.3806 

0,3839 

0.39378 

1 .6 

0.5 

0,5047 

0.51 676 

2.0 

0.6121 

0,6166 

0.62977 

2.4 

0,711 

0,7163 

0.72899 

2.8 

0,796 

0,8 

0.81152 

3.2 

0,864 

0,868 

0.87609 

3.6 

0,915 

0.919 

0.92333 

4.0 

0.951 

0,954 

0.95552 

4 .4 

0,975 

0.976 

0,97587 

4.8 

0.988 

0.990 

0,98779 

5.2 

0,989 

0,992 

0.994^ 

5.6 

0,993 

0,993 

0,99748 


aiAPTER 6 


CCBXUSIOMS AMD ^JGGEStlOflS 

1 ) Th® applicability of tho adaptaiicsi 

teehnlc^o ha® haen demon stra tod by the thxe© physical 
problem® ©onsidorod. the problems are two-Klimensimal 
but correapotid to diverse applications whidi proaaises the 
use of this technique to a vdde rarjgo of two-<ii»en»ioRal 
probletaa* It is easily possible to extend the method to 
three dimenslcfial cases also.. 

2) This technique la particularly useful in the 
case of non-liiioar problems 11 fee th<» test problem of 
chemically reacting lar4nar flows. Adaptaticsn is carried 
out after each iteratic4J to swcescively improve up«»i 
the prediction of solutitfi variables in the large 
gradient regl«»\s# 

3) The weight ccmstint shcsuld be a large value to 
retain the exiremum pr<^rtles of the elliptic generatimn 
system. 

If ^e extrtBum principles are violated the resultiiig 
grid may b« to© swBSh shewed and in some cases the interior 
nodes sheet out side the boundary wlii^ ie not ailcra»i* 



4) Fintlly it is Bhmm that this adaptive finite 

element teehniqee has several advantages over th® other 
tdapt3.ve grid methods as explained oel<*rt 

(a) The teohnicitj® is best suited for coupling 
vdth the finite element method uhish is used 
vddely novadays. Standard finite elawnt routines 
are used in the present i?ork and so tte software 
is very simple, 

(b) adaptation is carried out ail over the 
domain in me stroke unlike the line by line 

or series of one-cHiaensioreil adaptations folio*«d 
in the attrastion repulsion ith«e or the spring 
ana logy# So the gradients are sensed all over 
ttm domain at the same tiw and accordingly ^ 
point a ate moved* 

(«) The presentlv available grid adaptaticn 
raeti^ods are primarily meant for ftnite difference 
soiutiaift and so c-nmolicdted traiwforaatlon# of 
fwetnii"^ e<patiois» of the problems are necessary to 
oarry out aompittatleirts In the erthogoiial dotsain* 
mesh compien traaisf©r»tioiis art net in 

the mm of finite eleiwt’^# 

Cd) the autenaatio mesh rtlifieiisiils sthemes avails hli ^ 
for filnito eleasfits inareti# the ooapitatioiial 
tiiit for nniWfrlcal siwuM'lieii* Another dtsadvanteft ! 
«dltli thi^ schemes it that tlit grid it not dyntniitli; 

iis^ inierattioii it ntoetsarii ' 



The pErescsit seJ^aie is dym anica 11 y adaptive and 
it daes not refine the by adding sfl^re eleraenta. 
So ccfliputationai time done not increase veryimsh* 

fjjasHsrioiia 

1) in the present wcackt the xlrst derivative is 
takdi as the ccmtrol fyncti«t ’'Ith the sign of second 
derivative for oonicoliing the dir'?ction of node sawement* 
More detailed study is necessary for the selectiott of appro- 
priate ccntrol fync*Ucns* 

2) A ewrvatore term a$ indi^'^ted by Thostpson Usl 
cm be taken along 'ith the first derivative# to avoid 
clustering of nodes near sharp cs’'nors* 

3) In fast the original s»sh is not eoatable for the 
physical problem the adaptation my coapomd the ctfficnltiea 
and sivart osciliations may occur in the aoluticn* So 
techniqvea mst be developed to provide a good initial mesh* 

4) The constant *W* should be based on i<»0 

grid properties iwch as aayiMii^miniiiHffl spaoliig bet^iOin 
nodes or an*iiao8|^ad.niaim> angle betwen grid lines* 



APPENHIX A 


For adipotling the weah la a partictilar i^yalcal 
probleiB the prograntne is stored ip three files# for 
copvonleaee , The first file is tho original proble® coda 
which fimis the solutias of the problem for coarse grid* 

The solulicw is checked for cmvGrg(mce* If the adaptaticii 
is to be corriod out then subroutine Ll'ilffiHD in the seccmd 
file is called* This file adapts tii." mesh and if the new 
mesh ia accoptablf* it tails the INTPOt routine which Is in 
the same file. The routine TrlTPOt interpolatts the values 
of th© Yariabl®(s) to the new grid points# This inteipolatior 
is tarried out ch theifebasis of the previous grid and the 
pravioMs variable wines# These interpolated valttes beca» 
a part of Urn iterative procedure for the Min probJUiffi# Aft®\ 
this ©oeaeation the routine subsUtatas the new grid coordinat, 
in place of the previous value© and returns the new grid to 
the main pri^lem and again the solutlai procedure is carried 
The third fll® c®nstains the general routines used is finite 
elements lit® SHAPES# qjACOB, OTFItlt FMfS etc# sliich c«i 
referred in Qol * The flow chart for th® psm*Awm »ntiMe«! 
above can be seen in Appendix Q# : 




Flow chart for adaptation 
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